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Abstract 

In practical life, we encounter many systems that cannot be studied directly, either due to their high cost or because 
some of these systems cannot be studied directly. Therefore, we resort to the simulation method, which depends 
on applying the study to systems similar to real ones and then projecting these results if they are suitable for the 
real system. The simulation process requires a good understanding of probability distributions and the methods 
used to transform random numbers that follow a regular distribution in the field [0,1] into random variables that 
follow them, so that we can achieve the greatest possible benefit from the simulation process and obtain more 
accurate and appropriate results for all conditions that arise. In previous research, we presented a neutrosophical 
vision of the process of generating random numbers that follow a regular distribution in the field [0, 1] and some 
techniques used to generate random variables, such as the inverse transformation technique that was used to 
generate random variables that follow a uniform distribution in the domain [a, b] and the exponential distribution, 
the rejection and acceptance technique, which was used to generate random variables that follow the beta 
distribution, and the mixed technique, which was used to generate random variables that follow the Poisson 
distribution. In this research, we present a neutrosophic study to generate neutrosophic random variables that 
follow the gamma distribution, a distribution that is frequently used in engineering applications. 


Keywords: Simulation, Generating Random Numbers, Neutrosophic Logic, Generating Neutrosophic Random Numbers, 
Rejection and Acceptance Technique. 


1 | Introduction 


Operations research is the applied side of mathematics, and since its inception, it has helped improve the 
performance of many systems that use its methods in their workflow. Among the methods of operations 
research is the simulation method. The importance of the simulation process in all branches of science comes 
from the great difficulty that we can face when studying the operation of any system. From real systems, due 
to the high cost or the fact that some systems cannot be studied directly, the simulation process depends on 
generating a series of random numbers subject to a regular probability distribution in the range [0,1], then 
converting these random numbers into random variables subject to the law of probability distribution 
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according to which the system to be simulated operates, the attention of researchers and scholars interested 
in the simulation method has focused on providing scientific methods and methods that help in obtaining 
the transformation of random numbets into random variables that follow the probability distributions most 
used in practical applications, and in the face of the great revolution brought about by neutrosophic logic in 
All fields of science, which are confirmed by the research presented in this field [1-17], are a continuation of 
the research we presented previously, the aim of which is to present a neutrosophical vision of Monte Carlo 
simulation [18-22]. In this research, we present a neutrosophic study to generate neutrosophic random 
vatiables that follow the gamma distribution, a distribution that is frequently used in engineering applications. 


2 | Discussion 


In this paper we have founded a new type of hybrid algebra, called Neutrosophic TwoFold Algebra and its 
corresponding Neutrosophic TwoFold Law, by combining two algebras: 


2.1 | The Principle Behind the Simulation Process 


It is the generation of a series of random numbers that follow a uniform distribution over the field [0,1], and 
then converting these random numbers into random variables that follow the distribution according to which 
the system to be simulated operates. There are many techniques used to generate these random numbers, 
such as the mean square, the product mean., the Fibonacci method and others, and to generate neutrosophic 
random numbers we presented a study in the research [18] and we used the mean square method and we 
reached the following results: 


e The mean square method for generating a series of random numbers is defined by the following 
relation: 


Riga = Mid[R?] ; i = 0,1,2,3, ... 


Where Mid symbolizes the middle four ranks of R?, and R; is chosen, ie., a fractional random 


number consisting of four ranks (called a seed) that does not contain a zero in any of its four ranks. 


e We convert these random numbers into neutrosophic random numbers and here we distinguish three 
states of the field [0,1] with a margin of indeterminacy, in the three states we have € € [0,n] and 
0<n<l 


First case: [0 + €, 1] indeterminacy at the minimum of the range. We substitute in the relation: 


R R; —€ 
iNn~T_-¢ 
Second case: [0,1 + €] indeterminacy at the upper limit of the range. We substitute in the relation: 


ee 
ee 


Third case: [0 + €,1 + €] Indeterminacy in the upper and lower limits of the field. We substitute 
in the relation: 


Rig = Ry € 
2.2 | Generating Neutrosophic Random Variables 


To generate neutrosophic random variables that follow the probability distribution according to which the 
system to be simulated operates, we distinguish three cases: 


First case: Neutrosophic random numbers and probability distribution in the classical form. 


Second case: Classical random numbers and the probability distribution is in the neutrosophic form. 
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Third case: Neutrosophic random numbers and probability distribution in neutrosophic form. 
2.3 | Generating Neutrosophic Random Variables Based Gamma Distribution 


The gamma distribution is one of the important distributions used in engineering applications. Given this 
importance and so that we can provide more accurate simulation results suitable for all conditions, we present 
in this research a neutrosophic study of the process of generating random variables that follow the gamma 
distribution, based on the classical and neutrosophic studies presented in the references [18- 25]. 


2.3.1 | Classic formula for gamma distribution [23, 24] 


The probability density function is given by the following classical formula: 


a 
fa) = ate ; x>0,a>0,8>0 


T(x) = yo he “dx 
0 


where x => 0,a > 0,8 > 0 are non-negative values, 

Special cases of gamma distribution: 
e As the value of @ increases, the gamma distribution approaches the normal distribution. 
e When a = 1, the gamma distribution matches the exponential distribution. 


2.3.2 | Neutrosophic Formula for Gamma Distribution [25] 


If we take a or B or both, a neutrosophic number, where the neutrosophic number is written in the following 
standard form, N = a + bI, where a and b are real or complex coefficients. a represents the definite part 
and bI represents the indefinite part (indeterminacy) of the number N, and it is possible that [Ay Az | ot 
{Ay, Ay } or...otherwise it is any set close to the real value a,we obtain the neutrosophic gamma distribution. 
Here we distinguish the following cases: 

1. We take @ a neutrosophic value, B classic value, then the probability density function is written in the 

following form: 
Ee 

(ay) 
x20,ay >0,B >0 


xtn-1e-B.x 


fu) = 


2. We take B a neutrosophic value, a classic value, then the probability density function is written in the 
following form: 


fy (x) = By_ x 4-1 e-Bu x 


x20,a>0,By >0 
Thus, we obtain the following special cases of the neutrosophic gamma distribution: 
e As the value of @ increases, the gamma distribution approaches the neutrosophic normal distribution. 
e When a = 1, the gamma distribution matches the neutrosophic exponential distribution. 


3. We take a and f are neutrosophic values, that is, we take them in the form. Then the probability density 
function is written in the following form: 
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By’ 


xtn-19-BN-X 
ray) 


fn(x) = 


x20,ay >0,By > 0 


Where By =P +y anday=a+o6 and y,6 is the indeterminacy of the mediators B, a and it can 
be any set close to the real values B, a. 


2.3.3 | Generating Neutrosophic Random Variables Based Gamma Distribution 


From the classical study presented in references [23,24] and the neutrosophic study presented in the 
references [18,21,22], the general formulation of the technique used to generate neutrosophic random 
variables following a gamma distribution is as follows: 


1. 


3 | 


We search for the largest value of the function fy(x) on the domain of its definition by calculating the 
first derivative and making it equal to zero. We obtain the following value: 


an — 
= mod 


xn = 
N 


We replace this value with the function fy (x). We obtain M, the largest value of the function, which 
corresponds to the mode. Then we follow the following steps: 


a. We generate two random numbers from the range [0,1], such as Ry , Ro. 
b. We convert these two numbers into neutrosophic random numbers Ry, , Ry. 


; 1 
c. We form a new variable from one of the two numbers, for example xy = oe then the 
N1 


transformation field of x becomes [0, [, which is identical to the distribution field fy (x). 


( : ) 
__ fn Rag 
d. Wecompare Ry2 with —o 


tua) 
M 


to the gamma distribution. 


1 ; ’ 
e. IfRy2 > , we consider the variable Xy = 2, t° be a neutrosophic random variable subject 
N1 


1 : 
f. IfR, >x= Ra We reject the numbers Ry 1 ,Ryz and return to step (a). 
N1 


Example 


We explain the above through an example that corresponds to the third case of neutrosophic random 


num 


Find 


bers, and the gamma function is a neuttosophic function. 


a random variable that follows the gamma distribution defined by the following probability density 


function: 


By’ 


xtn-19-Bn-x 
r(@y) 


fu(x) = 


Where x >0,By € [0.2.0.25], ay € {3,5} 


Based on the data in the example, the study is divided into two parts: 


Dek 


Data 


| First Part 


x >0,ay =3,By € [0.02,0.025] 


First, we calculate M from the following relation: 
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ay —1 
By 


a las aoe € [8,10] 
Ny (02025| [o2zozs] ° 


My = 


a. We generate two random numbers that follow a uniform distribution in the range [0,1] using the mean 
square method given by the following relation: 


Ris1 = Mid[R?| ;i = 0,1,2,3, ... 
We take the seed Rg = 0.3176 and from it we get R; = 0.0869 ,Rz = 0.7551. 


b. We convert the two numbers into two neutrosophic random numbers on the range [0 + ¢€,1 + €] using 
the relation: 


Riy = Ri -€ 
We take in this example € € [0,0.03], we get the following two neutrosophic random numbers: 
Ry = 0.0869 — [0,0.03] € [0.0569,0.0869] 
Ray = 0.7551 — [0,0.03] € [0.7251,0.7551] 


c. We form a new variable from one of the two numbers, for example 


1 1 
= —_________€ [11,5075,17.5746] 


*N Rx, [0.0569,0.0869] 


fx{—— 
d. Tocompare Roy € [0.7251,0.7551] with inl), 
N 


ay 
e We calculate fy (—) we find: f(x) = Eu xen de But 


([0.2,0.25])3 
r(3) 
€ [0.1766,0.8031] 


funy) 


N 
e We calculate 
Mn 


fy ([11,5075,17.5746]) = ([11.5075,17.5746])2e~[0.2,0.25].[11.5075,17.5746] 


we find: 


pl 


fn (ene ) [0.1766,0.8031] 
y .1766,0. 
a (000 221;0.0803 
My [8,10] ' 


f a 
We have, Roy € [0.7251,0.7551] and ula) € [0.0221,0.0803] 
N 


We note that: [0.7251,0.7551] > [0.0221,0.0803] => Roy > 


Check the acceptance criterion, so we consider Xy = a € [11,5075,17.5746] 
N14 


To be a neutrosophic random variable subject to the gamma distribution. 


3.2 | Second Part 


Based on the data x >0,a@y =5,By € [0.02,0.025] 


We calculate M from the following relation: 
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ay —1 
By 
ay —1 5-1 4 


My = _ 5-1 4 
v= —B— = 02,025] ~ [02,025] © 11°70! 


Mn = 


a. We generate two random numbers that follow a uniform distribution in the range [0,1] using the mean 


square method given by the following relation: 
Ri41 = Mid[R?] ; i = 0,1,2,3, ... 
We take the seed Rp = 0.3176 and from it we get Ry = 0.0869 , Rp = 0.7551. 


b. We convert the two numbers into two neutrosophic random numbers on the range [0 + €,1 + €] using 


the relation: 
Rin = Ri -— € 
We take in this example € € [0,0.03], we get the following two neutrosophic random numbers: 
Riv = 0.0869 — [0,0.03] € [0.0569,0.0869] 
Roy = 0.7551 — [0,0.03] € [0.7251,0.7551] 


c. We form a new variable from one of the two numbers, for example: 


1 1 


= —_________¢ [11,5075,17.5746 
Ry1._ [0.0569,0.0869] ete; 


xn = 


fn 1 
d. Tocompare Roy € [0.7251,0.7551] with ——*+ Cr) 
Mn 


e We calculate fy (— -) we find: f(x) = Eu c@n-ter But 


([0.2,0.25])> 


l'(5) 
E [0.1766,0.8031] 


fy ({11,5075,17.5746]) = ([11.5075,17.5746])*e~[0.2,0.25].[11.5075,17.5746] 


fn (Cas ) tng ) [0.1766,0.8031] 
e We calculate ——“ we find: ——“++ = —_—_— € [0.011,0.04.01] 
My My [16,20] 


1 
We have, Roy € [0.7251,0.7551] and f(a) E [0.011,0.0401] 
Mn 


We note that: [0.7251,0.7551] > [0.011,0.0401] > Roy > 


Check the acceptance criterion, so we consider Xy = a € [11,5075,17.5746] 
N14 


To be a neutrosophic random variable subject to the gamma distribution. From the results of the study in the 


first and second sections, we note that: %y = —— € [11,5075,17.5746] 
N1 


Can be considered a neutrosophic random variable that follows a gamma distribution defined by the following 
probability density function: 


f(x) = Pu ytw-te-Pw* wh > 0,ay € {3,5} € [0.2,0.25] 
N = Ta ,whetre x 20,ay »t,Bn .2,0.25]. 
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3.3 | Classic Wording of the Example 


Find a random variable that follows the gamma distribution defined by the following probability density 


function: 

BF ied 

= —— ee —Bx 
FO) = Fayette 
Where x >0,a=3,8 = 0.2 
Datax >0,a=3,8 = 0.2 
First, we calculate M from the following relation: 
a-1 


B 0.2 O02 


a. We generate two random numbers that follow a uniform distribution in the range [0,1] using the mean 


square method given by the following relation: 


We take the seed Rp = 0.3176 and from it we get Ry = 0.0869 ,Rz = 0.7551. 


b. We form a new variable from one of the two numbers, for example R1, we find: 


Ri41 = Mid[R?] ; i = 0,1,2,3, ... 


1 


=—= = 11.5075 
~=R, 0.0869 


fle) 


c. Tocompare Rp with— 


— - — Be a-1,-Bx 
We calculate f (= ) we find: f(x) = T@* & 
(0.2)3 ; aes 
fy(11.5075) = ra) (11.5075)2e~92-115075 — 0.1766 

f= f= 

We calculate fa) we find: a) — 22766 — 0.0221 
M M 8 
1 

_ se) | 

We have, R2 = 0.755 land —— = 0.0221 


f 
We note that: 0.7551 > 0.0221 => R, > —** 


ee a 
Check the acceptance criterion, so we consider: Xy = - = 11.5075 
1 


To be a neutrosophic random variable subject to the gamma distribution. 


From the results of the study in the first and second sections, we note that: 


1 
ky = p= 11.5075 
1 


Can be considered a neutrosophic random variable that follows a gamma distribution defined by the 


following probability density function: 
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pr 
— " ya-1,-Bx 
f@ = par" 
Where x > 0,a = 3,8 = 0.2. 


The classic example is treated in the same way as before if the data are: 


Find a random variable that follows the gamma distribution defined by the following probability density 
function: 


pate ; x>0,a>0,8 > 0, wherex >0,a=5,8 = 0.2. 
T(@) 


4 | Conclusion and Results 


In this research, we presented a neutrosophic study, the purpose of which is to obtain neutrosophic random 
variables that follow the gamma distribution, a distribution that has many uses in engineering applications. 
We found through this study that the neutrosophic values of the required random variables have a margin of 
freedom and take into account all the conditions that the system’s operating environment passes through. 
What is intended to be simulated works according to the gamma distribution, while the classical study gives 
a specific value for the random variable that we need for the simulation process, a value appropriate to 
conditions similar to the conditions in which such data were adopted, and the results that we obtained through 
the neutrosophic presentation of the example and the classical presentation confirm these results, and 
therefore to obtain more accurate simulation results that take into account all conditions, we recommend 
using the neutrosophic study. 
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